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ABSTRACT 


On  the  basis  of  the  metric  in  a  gravitational  field,  a  two-way  Doppler 
shift  is  obtained  for  an  arrangement  in  which  a  signal  sent  from  a  ground- 
based  transmitter  is  received  by  a  satellite,  retransmitted,  and  picked  up 
by  a  receiver.  The  geometry  used  is  quite  general,  and  the  formula  obtained 
is  correct  to  order  (v/c)°  included.  When  both  transmitter  and  receiver  are 
on  the  earth,  much  of  the  general  relativistic  effect  is  seen  to  drop  out. 
Various  limits  of  the  Doppler  shift  formula  are  considered. 
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I.  INTRODUCTION 


This  report  derives  a  relativistic  Doppler  shift  formula  under  the  most 
general  condition  possible.  It  is  assumed  that  both  transmitter  and  receiver 
are  moving  in  an  arbitrary  direction  with  rates  not  necessarily  constant.  The 
derivation  is  done  with  the  use  of  general  relativity  and  includes  the  correc¬ 
tions  due  to  the  gravitational  field  on  the  motion  of  light  waves.  The  situation 
envisioned  is  one  in  which  a  transmitter  sends  a  monochromatic  light  signal 
to  a  receiver.  The  receiver  picks  up  the  signal,  multiplies  it  by  a  fixed  ratio, 
and  sends  it  back  to  a  second  receiver.  This  is  a  two-way  Doppler  arrange¬ 
ment.  The  second  receiver  may  or  may  not  oe  at  the  same  location  as  the 
transmitter;  at  any  rate,  the  geometry  is  left  to  be  completely  general. 

The  formula  derived  here  can  be  used  to  calculate  the  range  rate  of  a 
satellite,  in  which  case  the  transmitter  is  a  ground -stationed  antenna,  the 
receiver  the  orbiting  satellite,  and  the  second  receiver  a  tracking  antenna. 

The  high  degree  of  precision  of  the  formula  obtained  in  this  report  may  be 
necessary  for  the  correct  tracking  of  satellites  orbiting  at  large  altitudes, 
in  particular  at  synchronous  altitudes. 

The  one-way  Doppler  shift  is  derived  first;  the  derivation  is  subsequently 
extended  to  the  two-way  arrangement. 


II.  DERIVATION  OF  ONE-WAY  DOPPLER  SHIFT 


A.  METRIC  IN  A  GRAVITATIONAL  FIELD  AND  PATH 
OF  A  PHOTON 

The  general  form  of  the  space -time  relativistic  metric  is 

ds2  =  dx^  dxv  (1) 

where  the  quantity  ds  is  called  the  proper  time. 

In  the  case  of  a  gravitational  field  it  a  point  ~r  outside  a  sphere  of  mass 
M  whose  center  is  at  r  =  0,  the  metric  is  given  by  the  Schwartschild  space - 
time.  Its  form  is 


ds 


2 


(i  -  ^>dt2 . 4 

c 


dr 

1  -  (2m/ r) 


+  r2de2  +  r2sin2e  d2<t> 


(2) 


where  the  constant  m  is  defined  as 


m 


G  Me/c2. 


(3) 


In  Eq.  (3),  G  is  the  gravitational  constant,  M  the  mass  of  the  earth,  and  c 

c 

the  velocity  of  light.  The  geodesics  of  the  electromagnetic  radiation  is  found 
by  differentiation  with  respect  to  some  non-zero  parameter  |Jt.  (For  light 
ds2  =  0. )  One  has1 


*G.  McVittie,  General  Relativity  and  Cosmology,  Univ.  of  Illinois 
Press,  Urbana  (1965),  p.  80 


The  following  relation  also  holds 


g 


vv 


2 


0 


(5) 


Only  four  of  the  relations  defined  by  Eqs.  (4)  and  (5)  are  independent. 

Taking  the  equations  that  are  obtained  from  setting  X.  =  1 , 3,  4  (i.  e. ,  xV  =  t,  8,  <|>) 
and  using  Eq.  (5),  we  obtain  for  the  motion  of  electromagnetic  waves 


We  are  interested  in  motion  in  a  plane  passing  through  the  origin  r  =  0  that 
traces  an  orbit  about  this  point,  and  we  obtain 

6  =  tt/2  ,  d6/d  \i  =  0.  (10) 

This  boundary  condition  is  obviously  a  solution  of  Eq.  (8).  The  other 
equations  reduce  to 


/4  2m.  dt  _ 

(1  — — )  =  P  »  a  constant 


(id 
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_i  (r2  M)  =  0 

d(i  '  d|i' 


r2p-  =  h 

d|i 


a  constant  (12) 


c 


Eqs.  (11)  and  (12)  yield 


_ £___  *4  =  h=h* 

1  -  (2m7r)  dt  “  p 


Setting  fi  =  t  and  substituting  Eq.  (14)  into  Eq.  (13)  gives 


dr  _  .a  2m»  2  h 

dt  r  '  r  r 


(1  -S2) 

r 


This  equation  represents  the  radial  motion  of  photons.  The  quantity  h  is  a 

❖ 

constant  along  any  path;  however,  the  value  of  h  is  different  for  different 
paths.  For  radial  propagation  it  is  clear  that  h  =0,  since  d<J>/dt  =  0. 

B.  ONE-WAY  DOPPLER  SHIFT  CON  FIGURATION 

A  transmitting  station  T  is  fixed  on  the  moving  earth.  The  space-time 
coordinates  describing  the  transmission  of  an  electromagnetic  wave  are  given 
by  the  pair  of  events  (rt,  tt/2,  4>t>  tt)  and(rt>  tt/2,  4>t  +  A4>t»  tt  +  At^).  The  change 
in  the  4>  coordinate  A<j>.  is  due  to  the  rotation  of  the  earth,  and  At  is  the  time 

*  L 

elapsed  during  the  emission  of  the  wave.  The  coordinates  of  the  reception 

of  this  wave  on  the  satellite  S  define  two  events  (r  ,  tt/2,  <j>  ,  t  )  and  (r  +  Ar  , 

s  s  s  s  s 

tt/2,  <j>  +  A<t>  ,  t  +  At  ).  The  position  of  the  satellite  changes  by  an  amount 

8  S  S  S 
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Ar  and  A<J>  in  the  time  At  ,  which  measures  the  change  of  the  coordinate 
s  s  s 

time  during  the  arrival  of  the  wave.  By  comparing  the  time  intervals  At^  and 
Atg,  we  can  calculate  the  range  rate  of  the  satellite. 

C.  INVARIANT  DEFINITION  OF  FREQUENCY  AND 
FREQUENCY  SHIFT 

The  frequency  f^  of  the  wave  emitted  by  the  ground  transmitter  can  be 
defined  in  an  invariant  way  by 


f. 

L 


(16) 


where  As^  is  the  invariant  proper  time  interval.  The  frequency  fg  =f^  +  Af^ 
of  the  wave  emitted  by  the  ground  and  received  by  the  satellite  can  be  similarly 
defined  in  an  invariant  fashion  by 


f  =  f  +  Af 
s  t  t 


(17) 


The  two  proper  times  As  and  As  are  defined  in  Section  D.  Using  Eqs.  (16) 

t  s 

and  (17),  we  can  write 


ft  +  *ft 


(18) 


The  frequency  shift  Af^  will  be  shown  to  be  related  to  the  range  rate  of  the  satellite. 
D.  PROPER  TIMES 


In  the  geometry  we  are  considering,  the  proper  time  ds  is  defined  by 


ds2  =  (1  .121)  dt2  -  -V 

r  C6 


l.  2m  \_1  j  2  2  ,,2 

^  -  —)  dr  +  r  d<(> 


(19) 
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We  can  write  this  as 


[1  -  (2m/r)1dt2 
ll  -  (2m/r)]dt^ 

or 


ds 


,  2  .  2  ,  .2 
dr  +  r  d<(> 


ds 


2 


J_  }[1  -  (2m/ r)]-1  dr^  +  r2 
c2  [1  -  (2m/r)]dt2 


If  we  define  a  velocity  relativistically  as 


q 


1 

(1  -  2m/r)1/2 


[i  -  (2m/r)]_1 


» 


ve  obtain 


Using  this  form  for  the  definition  of  the  proper  time,  we  can  write  for  the 
definition  of  the  proper  times  As^  and  Asg 


2 


We  are  keeping  d0/dt  =  0  in  this  discussion. 


and 


As 


s 


(25) 


E.  COMPARISON  OF  THE  COORDINATE  TIME 

INTERVALS  At  AND  At 
_ t _ s_ 

To  be  able  to  compute  the  ratio  f  /£,.  of  Eq.  (18),  we  need  to  develop  an 

r  t 

expression  for  Af  /At  .  This  task  can  be  simplified  by  reference  to  Figure  i, 
c  s 

which  relates  the  coordinates  describing  the  emission  and  reception  of  the 
electromagnetic  waves. 

The  start  and  the  end  of  a  wave  are  signaled  by  photons  emitted  at 

and  r  +  Ar  ,  respectively.  These  photons  travel  radially  to  r  and  r  +  Ar  , 
t  t  s  s  s 

respectively.  For  radial  propagation  we  have 


If  =  c(l  -  2m/ r) 


(26) 


Then  for  the  first  photon  we  have 


dr 

1  -  2m/ r 


(2?) 


For  the  photon  signaling  the  end  of  the  wave  we  can  write 


<ts  +  A  V  -  (tt +  =  * f 


r  +Ar 
s  s 


dr 

1  -  2m/ r 


(28) 


rt+Art 


Subtracting  Eq.  (27)  from  Eq.  (28),  we  obtain 


At  -  At 
s  t 


■*/ 


r  +  Ar 
/•s  s 


dr 

1  -  2m/ r 


rt+Art 


r 

.  fs _ 

J  1  -  2m/  r 

rt 


With  reference  to  Figure  1  we  can  rewrite  Eq.  (29)  in  the  form 


At  -  At 
s 


-  if  (  dr  f 

t  c  J  1  -  Zm/r  J 
.  path  1  |r 


r  +Ar  I 
1  s  s' 


dr 

1  -  2m/ r 


+Ar  -6r  i 
s  s  s> 


f  1  dr  /  dr 

"  J  1  -  2m/ r  ~  /  1  -  2m/ r 

|r  |  path  2 

L 


In  the  evaluation  of  the  integrals  in  Eq.  (30),  it  is  useful  to  bear  in  mind  that 

Ar  and  Ar  are  the  distances  travelled  by  transmitter  and  satellite  during  the 
t  s 

time  of  transmission  and  reception  of  one  cycle  of  the  signal.  Using  f  a  1.8 

Q  _7  _A  t 

X  10  cps  we  get  Ar  »  2.  6  X  10  m  and  Ar  »  1.7  X  10  m  for  satellites  at 

t  s 

synchronous  altitudes.  This  means  that  paths  1  and  2  are  parallel  to  a  very 
high  degree  of  accuracy.  Furthermore,  the  greatest  error  we  can  make  in 
saying  that  the  lengths  of  paths  1  and  2  are  the  same  can  be  found  by  letting 
Ar^  =  0  and  Ar^  a  2  X  10  ^  m  and  assuming  that  Arg  is  normal  to  path  2.  With 
rj  and  r^  designating  the  lengths  of  paths  1  and  2,  respectively,  we  have 


/  2  ,  A  2  ,  1  Ar 

I  r2  +  Ar  a  r2  +  j  — 


Ar  «  r. 
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-19 

With  the  satellite  at  synchronous  altitude  we  obtain  <  10  m 

remaining  integrals  in  Eq.  (30)  are  of  the  order  of  Art  and  Arg;  thus, 
the  integrals  over  paths  1  and  2  equal  amounts  to  rejecting  quantities 

3 

Oi.Ar  ),  which  is  perfectly  justified.  We  now  have 


At  -  At 
s  t 


Uf  . 

dr  r  _ 

dr  \ 

c\  J 

1  -  2m/r  J  1  -  2m/r  I 

Vbr 

6r. 

/ 

This  gives 


At  -  At 
s  t 


1  /  6rs 

c  \  1  -  2m/ r 

\  s 


1 


6rt  \ 

-  2mTrt  J 


With  reference  to  Figure  1  we  can  write 


6r  =  Ar  cos  6  =  Ar 

s  s  s  s 


and 


6r 


t 


Art  cos  et  =  Art 


Substituting  Eqs.  (34)  and  (35)  into  Eq.  (33),  we  obtain 

Ats  -  Att  =  c  (^'A?s  -  S*A?t) 


.  The 
setting 
of  order 


(32) 


(33) 


(34) 


(35) 


(36) 
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where  we  have  let 


A  = 


1 

1  -  2m/ r 


r  -  r, 
s  t 


S  |7-Tt| 


(37) 


and 


B  1  -  2m/ r. 


r  -  r, 
s  t 


*  'rs  -  rtl 


(38) 


Eq.  (36)  can  also  be  written  as 


i  /_►  Ar  Ar  At  \ 

At  -At  =  ±  A-— ^  At 

s  t  c  \  At  At  At.  /  s 

'  s  s  t  / 


(39) 


or 


i  .  Att 

Ats  "  Att  =  c  Avs  -  B' 


■ ,  Att  \ 

l-vt  A rj 


At 


(40) 


Dividing  by  At  gives 
s 


i  - 


Att  i  /-  _  -  ^  Att  \ 
xri-  =  -  A  v  -  B-v.  ~) 
At  c  \  3  t  At  / 


(41) 


Rearranging  terms  we  obtain 


Att  B  vt  Att 
At  "  c  ‘  At 


A*  v 


=  1  - 


(42) 
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or 


At 

s 


1  -(A- "v*  /c) 
s 

1  -(B-vt/c) 


Writing  A  and  B  explicitly,  we  have 


(43) 


which  is  the  relation  between  At  and  At  we  sought. 

t  s 

F.  ONE-WAY  DOPPLER  SHIFT  FORMULA 

We  can  now  gather  Eqs.  (19),  (24),  (25),  and  (44)  to  obtain 


(45) 


This  expression  is  relativistically  correct  to  the  order  of  the  terms  neglected 

in  the  calculation  of  Eq.  (44).  The  terms  neglected  are  of  order 
vs  2  v  4 

—  (Ar)^  =  (— )  .  To  the  same  accuracy  we  can  replace  qt  with  vt  and  qg  with 

v  .  Thus  we  have 
s 


-13- 


in.  TWO-WAY  DOPPLER  SHIFT  FORMULA 


To  obtain  the  two-way  Doppler  formula,  we  need  a  relation  between  the 

frequency  transmitted  by  the  satellite  ftg  and  the  frequency  received  by  the 

ground-based  receiver  fr-  If  the  satellite  transmitter  sends  out  a  signal  with 

a  frequency  k  times  the  received  frequency  f  ,  we  have  f.  =  kf  ,  so  that 

s  t  s  s 


f 

r 


f 

s 


f 

s 


(47) 


No  new  derivation  is  needed  for  the  calculation  of  f  /L  .  We  need  only 

r  ts  7 

observe  that  the  analysis  carried  out  before  still  applies  if  we  replace  the 
subscript  t  with  s  and  the  subscript  s  with  r  in  the  righthand  side  of 
Eq.  (46).  We  now  have 


(48) 


The  final  result  is  now  obtained  by  substitution  of  Eqs.  (46)  and  (48)  into 
Eq.  (47): 


(Con't) 
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WJ  ^IIWP 


This  formula  is  relativistically  correct  to  ojjv/c^ 

We  will  assume  that  the  transmitter  and  receiving  stations  are  ground 
based  and  located  at  the  same  geographical  coordinates.  Then  we  can  set 


r  =  r  and  v  =  vA 
r  t  r  t 


(50) 


but 


?r  f  rj.  and  Vj.  t  vj.  .  (51) 

The  vectors  defining  the  position  and  velocity  of  the  receiver  and  transmitter 
are  not  the  same,  because  during  the  time  interval  necessary  for  a  photon  to 
travel  from  the  earth  to  the  satellite  and  back,  the  earth  has  rotated  and  the 
direction  of  the  vectors  ~r  and  *v  has  changed.  The  extent  to  which  Eq.  (49) 
is  influenced  by  this  effect  will  be  calculated  later. 

Using  Eq.  (50)  we  can  simplify  Eq.  (49)  to  yield 


-16- 


where  we  have  defined 


It  should  be  noted  that  in  the  two-way  Doppler  effect,  in  which  transmitter 
and  receiver  are  both  on  the  earth,  the  general  relativity  effects  of  the  order 
of  2m/;  have  dropped  ^ut  and  they  enter  only  as  (2m/r)(v/c),  which  is  one 
order  o 2  magnitude  smaller  than  2m/ r. 

It  is  useful  to  know  the  order  of  magnitude  of  the  quantities  in  the  right- 
hand  side  of  Eq.  (52).  With  the  assumption  that  the  satellite  is  at  geosynchron¬ 
ous  altitude,  we  have 


1.5  X  10"6 


2m 

r 

s 


2  X  10 


-10 


2m  „ 


2S  1.4  X  10* 


(54) 


Therefore,  if  we  are  to  expand  Eq.  (52)  to  terms  of  order  (v/c)3,  we  must 
retain  the  terms  of  the  form  (2m/r)(v/c). 

We  noted  previously  that  Eq.  (52)  is  valid  only  to  terms  of  order  of 
(v/c)^  included.  Treating  the  scalar  products  in  Eq.  (52)  as  small  in  compari¬ 
son  to  1,  expanding  in  powers  of  v/c,  and  retaining  only  quantities  of  the  order 
of  (v/c)3  gives  the  following  expression. 


k  1  +-Mv-R  +T.R  -■/  •  R  -  v*  •  R.  \  + 
c\ss  tt  rs  st/ 


♦  7  [( v  Kf  *  (v  *,)2  *  (v  MfV  *.)  -  (v  *.)  •  (y  *,)(V  K) 

■K- K)  'ft'  N)(V  *t)*(v  KK- *.)]♦  (Con-,) 
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The  limit  of  this  equation  when  we  set  rf  =  rt  and  =  vt  is  examined  in  the 
the  Appendix. 

We  can  obtain  a  better  understanding  of  the  contribution  of  each  term 
on  the  righthand  side  of  Eq.  (55)  by  looking  at  the  relationships  between  the 
quantities  r*,  v^,  and  vr.  We  begin  by  noting  that  for  a  point  on  the  earth, 

moving  with  angular  velocity  co. 


V  =  co  x  r 


(56) 


and 


r(t  +  At)  »  r(t)  +  vAt  =  ?(t)  +  co  X  r(t)At  .  (57) 

Therefore  we  can  write 


yj.  =  w  X  rt  (58) 

v  =  w  x  T  (59) 

r  r 
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rt  +  u>rtAtvt 


(60) 


r  =  r.  +  ^  X  r.  At  = 
r  t  t 


It  is  now  possible  to  obtain  a  relationship  between  the  unit  vectors  Rg  and  R 

By  definition  of  R  and  using  Eq.  (60),  we  have 
s 


R 

s 


?  -  ?  rx  +  oar  A  tv  -  ? 

r  s  t  t  t  s 


I  ?  -  ?  I  I  ?  +  oar  Atv,  -  r*  | 
r  s  t  t  t  s 


(61) 


The  relative  importance  of  the  second  term  on  the  righthand  side  of  Eq.  (61) 

is  |u>rt  At  ^ l/f^t  I  =  -  3  X  10”^,  of  the  order  of  v/c.  The  value  of  At 

used  to  estimate  oaAt  is  the  time  necessary  for  a  light  wave  to  go  from  the  earth 

to  the  satellite  at  synchronous  altitude  and  back,  oa  is  the  angular  velocity  cf 

the  earth.  Since  oar  At  is  small  compared  to  r  and  r  ,  a  fruitful  way  of 

t  t  s 

writing  Rg  is 


r.  -  r 
t  s 


r  +  oar.  Atv, 
s  t  t 


oar^Atv^ 
r  +  oar  AtvJ 


(62) 


The  absolute  value  in  the  denominator  c  m  be  expanded  as  follows: 


|rl  -  T  + 
1  t  s 


r  +  oar  AtvJ 
s  t  t 
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IV.  CONCLUSIONS 


The  Doppler  formula  obtained  in  this  report  is  a  general  relativistic 
two-way  Doppler  shift  for  a  signal  emitted  by  a  ground  transmitter  and 
received  by  a  ground  receiver  after  being  reflected  by  an  orbiting  satellite. 
The  analysis  is  carried  out  under  most  general  conditions,  but  the  approxi¬ 
mation  made  in  the  computation  of  the  integrals  needed  to  evaluate  the  ratio 

At  /At  makes  the  result  valid  to  quantities  of  order  (v/c)^  when  the  satellite 
t  s 

is  at  synchronous  altitude.  The  formula  obtained  is  sufficiently  accurate  to 
evaluate  satellite  range  rate  to  (v/c)'5  c  =  10  cm/sec. 
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APPENDIX 


In  this  Appendix  we  want  to  analyze  the  limiting  value  that  Eq.  (55) 

assumes  when  we  set  ?  =  r.  and  v  -  v.,  From  the  definitions  of  R,  and  R  , 

r  t  r  t  t  s 

it  is  clear  that  setting  r^  =  r^  implies  R^  =  -Rg  =  R.  R  is  a  unit  vector  along 
the  direction  joining  the  earth-based  stations  and  the  satellite. 

A  little  algebra  shows  that  the  terms  in  c“*  become  -2(v  -vJ*R>  the  terms 
in  c“2  become  2[(vg  -  v^)*  R]  ,  and  the  terms  in  c"-3  reduce  to  -2[(vg  -  v^)*  R] 

+  2(v*R)^(v  •  R)  -2(v  •  R)  -  (v  •  R)  (4m/r  )c2  +  (v  •  R)  (4m/r  )  c^.  Using  these 

L  S  S  S  St  t 

results  in  Eq.  (55)  gives 


f 

r 


(  2(v  -  v  )•  R  2[(vs  -  vt).  R]2  2((v  -  vt)-  R]3 

kJ* - - -  +  - 2 - 3 -  + 

(  c  c 

+  4  |^2(vt-  R)2(vs-  R)  -2(vg.  R)2(vt-  R)  -(v.- R)^.  c2  +(vt-  R)*j“CZj  j 


(A-l) 


Defining  vg-  R  =  ug  and  v  •  R  =  we  can  rewrite  Eq.  (A-l)  in  the  following 
form 


f 


r 


=  k 


2(us  -  ut)  2(us  -  ut)  2(us  -  utr 

- - -  +  —  "  - . A -  -  - 5 - 


9  2 

2u  u 
t  s 


9  2 
2u  u. 
s  t 

— 3— 


4m  us 

r  c 
s 


4m  ut 


(A-2) 


We  now  define  the  range  rate  r  to  be  the  rate  of  the  separation  between  the 
transmitter  frame  and  satellite  frame,  as  seen  by  an  observer  on  the  trans¬ 
mitter  frame.  Because  of  the  law  of  addition  of  velocity  in  relativity,  it 

should  be  clear  that  r  *  u  -  u. ;  the  correct  relation  is,  in  fact, 

s  t 
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r  = 


u  -  u, 
s  t 


1  -  <UsUt/c  > 


(A-3) 


We  can  rewrite  this  expression  more  conveniently  for  our  purpose  as 


u  -  u 
s  t 


1  UsUt\  • 

1  )r 

k  C  / 


(A-4) 


or 


u  u, 

•  s  t  • 

u  -  u^  =  - - =-  r 

st  2 

c 


(A-5) 


To  orders  of  (u/c)3  we  have 


us  _  Ut  =  (Us  "  Ut)  +  °[<u/c>4] 


(A-6) 


or 


2  2 

u  uA  u,  u 
•  s  t  .  t  s  ,  ~ 

u  -  u .  =  r - 5“  + - y-  +  O 

s  t  2  2  | 

c  c  L 


j^(u/c)4j 


(A-7) 


Substituting  Eq.  (A-7)  into  Eq.  (A-2)  and  keeping  only  terms  of  the  order  of 
(u/c)3,  we  obtain 


9  2  9  2 

2u  u,  2u^  u 


•  <:u  u.  cu.  u  .  •  .  £  ,  •  x3 

T2*^)  -<f)  * 


,  2  ,  2 

.  2utus  2usut  4m  us  ,  4m  ut 

+  — 2 - F~  -T"  T  +  T_  T 

CCS  t 
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or 


f 

r 


(A-9) 


The  standard  formula  derived  for  a  two-^ay  doppler  shift,  which  is  valid 
for  near-earth  geometries  and  does  not  include  general  relativity  effects,  is 


(A  -10) 


where  the  notation  (  )'  is  used  to  differentiate  from  Eq.  (A-9). 
By  expanding  Eq.  (A -10)  to  orders  of  (r/c)^,  we  obtain 


(A-l  1) 


Comparison  of  Eqs.  (A-9)  and  (A —1 1 )  shows  that  they  agree  to  order  (r/c)^, 
the  only  difference  being  the  terms  and  -  ^£1  ~ ,  which  depend 

exclusively  on  general  relativity  effects. 
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II.  DERIVATION  OF  ONE-WAY  DOPPLER  SHIFT 


A.  METRIC  IN  A  GRAVITATIONAL  FIELD  AND  PATH 
OF  A  PHOTON 

The  general  form  of  the  space-time  relativistic  metric  is 

d®2  =  8^  dx**  dx*'  (1) 

where  the  quantity  ds  is  called  the  proper  time. 

In  the  case  of  a  gravitational  field  at  a  point  T  outside  a  sphere  of  mass 
M  whose  center  is  at  r  =  0,  the  metric  is  given  by  the  Schwartschild  space - 
time.  Its  form  is 


ds2  =  (1  -  ^p)dt2  -  -J 

c 


dr 

1  -  (2m/r) 


+  r2d62  +  r2sin20  d2«t> 


(2) 


where  the  constant  m  is  defined  as 


m  =  GMe/c2.  (3) 

In  Eq.  (3),  G  is  the  gravitational  constant,  M  the  masg  of  the  earth,  and  c 
the  velocity  of  light.  The  geodesics  of  the  electromagnetic  radiation  is  found 

by  differentiation  with  respect  to  some  non-zero  parameter  p.  (For  light 

2  2 
ds  =  0.)  Following  McVittie  and  Schusterman  one  obtains  the  following 

expression; 


*G.  McVittie,  General  Relativity  and  Cosmology,  Univ.  of  Illinois 
Press,  Urbana  (1965),  p.  &0 

2L.  Schusterman,  Relativistic  Red  Shift  and  Its  Application  to  the  Ranee  Rate 
Measurements  of  Satellites.  LMSC/584748,  Lockheed  Missiles  and  Space  Co,  , 
Sunnyvale,  Calif.  (11  May  1967). 
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The  following  relation  also  holds 


4 


2 


g 


vv 


2 


0 


(5) 


Only  four  of  the  relations  defined  by  Eqs.  (4)  and  (5)  are  independent. 

Taking  the  equations  that  are  obtained  from  setting  X  =  1 , 3,  4  (i.  e. ,  xV  =  t,  0,  <|>) 
and  using  Eq.  (5),  we  obtain  for  the  motion  of  electromagnetic  waves 


We  are  interested  in  motion  in  a  plane  passing  through  the  origin  r  =  0  that 
traces  an  orbit  about  this  point,  and  we  obtain 

0  =  tt/2  ,  dG/dp  =  0.  (10) 

This  boundary  condition  is  obviously  a  solution  of  Eq.  (8).  The  other 
equations  reduce  to 


(1  -  T1)  =  p 


a  constant 


(id 
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